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Abstract 

In this paper we prove an approximate formula expressed in terms of elementary functions for the 
implied volatility in the Heston model. The formula consists of the constant and first order terms 
in the large maturity expansion of the implied volatility function. The proof is based on saddlepoint 
methods and classical properties of holomorphic functions. 

1 Introduction 

In financial markets stochastic models are widely used by traders and risk managers to price and hedge 
financial products. The models are chosen on economic grounds to reflect the observed characteristics 
of market data, such as leptokurtic returns of asset prices and random instantaneous volatility of these 
returns. Pricing and hedging in realistic models are numerically intensive procedures that need to be per- 
formed very quickly because future investment decisions depend on the outcomes of these computations. 
Therefore market practitioners pay particular attention to tractability when deciding which model to use. 
Among the plethora of possible choices, stochastic volatility models are extremely popular particularly 
in equity, foreign exchange and interest rate markets since (i) they feature most of the characteristics of 
these markets, and (ii) they are numerically tractable. 

In practice, stochastic volatility models are first calibrated on market data, then used for pricing. 
Pricing financial products is mathematically tantamount either to solving a PDE problem with boundary 
conditions (the final payoff of the product) or to calculating the expectation of this payoff using prob- 
abilistic tools such as Monte Carlo simulation or stochastic integration. For most models closed-form 
formulae are scarcely available, and accurate algorithms have been extended and used such as finite- 
differences [23], ADI schemes [3T], or quadrature [5] methods. The calibration step involves a proper 
selection of the data to be fitted by a model. A common practice is to calibrate the so-called implied 
volatility rather than option prices directly. The implied volatility is a standardised measure of option 
prices which makes them comparable even though the underlying assets are not the same. Since this 
calibration step is based on optimisation algorithms, the lack of a closed-form formula for the implied 
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volatility makes it very time consuming. For instance, the SABR stochastic volatility model has become 
very popular because a closed-form approximation formula for the implied volatility was derived in |18j 
and hence made the model easily tractable. Likewise, perturbation methods as developed in |13j have 
proved to be very useful for obtaining a closed-form approximation formula of option prices. Although 
these methods only hold under some constraints on the parameters, they provide useful initial reference 
points for calibration. 

The Hcston model [5D] introduced in 1993 has become one of the most widely used stochastic volatility 
models in the derivatives market (see P3], [37], [2J, [3], [IB])- In this paper, we provide a closed- form 
approximation for the implied volatility in this model. The idea behind this result is the following: 
suppose one wants to calibrate the Hcston model on market data. This can be performed in two different 
ways: (i) one uses a global optimisation algorithm, which involves computing the implied volatility at 
each observed point until the algorithm converges; (ii) one specifies an initial set of parameters for the 
model and runs a local optimisation algorithm such as the least-squares method. The latter solution is 
the most widely used in practice since it is less computer-intensive. However its robustness heavily relies 
on the initial set of parameters to be specified. Simple closed-form approximations for the Heston model 
make this choice robust and accurate. One first calibrates the approximation on market data — which is 
straightforward since this is a closed-form — then one uses this calibrated set of parameters as a starting 
point in the whole calibration process. 

Let us consider an European option with maturity t and maturity-dependent strike K = So exp(xt), 
then our main result is the following asymptotic closed- form formula for the implied volatility a\ (x) : 

as the maturity t tends to infinity, where cr^ is defined in (TT8"]) . A in (|12p and Abs in (TT71) . For a constant 
strike K = Sq exp(x), we obtain the following formula: 

ai{x) = 8V*(0) + t~H (x(2p*(Q) - 1) - 2 log (-A(0)V2V*(0))) +o(t~ 1 ) (2) 

as the maturity t tends to infinity, where V* is given by (JHJ) and p* by ([9]). 

It is a well-known fact that for a fixed strike, the implied volatility flattens as the maturity in- 
creases [31] : this is confirmed by formula ([2J above, the zeroth order term of which was already known 
(see [27], [12]). However, the maturity-dependent strike formulation in formula ([T]) above reveals that 
the implied volatility smile does not flatten but rather spreads out in a very specific way as the maturity 
increases. 

In the fixed-strike case, Lewis [37J pioneered the research on large-time asymptotics of implied volatility 
in stochastic volatility models by studying the first eigenvalue and eigenfunction of the generator of the 
underlying stochastic process. Recently Tehranchi |33j studied the large-time behaviour of the implied 
volatility when the stock price is a non-negative local martingale and obtained an analogue of formula ([2]) 
in that setting. Comparatively, there has been a profusion of work on small-time asymptotics, based 
on differential geometry techniques [19], PDE methods [6] or large deviations techniques ([11] and [TO]). 
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Likewise, many papers have studied the behaviour of the implied volatility smile in the wings (see [I], 

0, m, EH, EU). 

The proof of our main result, Theorem 12.21 is based on two methods: first, we use saddlepoint 
approximation methods to study the behaviour of the call price function as an inverse Fourier transform. 
This idea has already been applied by several authors, including [pj, fJ5], [T] and [30] m order to speed 
up the computation of option pricing algorithms based on inverse Fourier transforms. We are also able 
to obtain the saddlepoint in closed form, thus avoiding any numerical approximations in determining it. 
The second step in our proof relies on Cauchy's integral theorem and contour integration for holomorphic 
functions in order to obtain precise estimates of call option prices in the large maturity limit. 

The paper is organised as follows. Section[2]contains the large-time asymptotic formula for call options 
under the Heston and the Black-Scholes models, both in the maturity-dependent and in the fixed-strike 
case. The proof of the main theorem, Theorem 12.21 is given in Section In section [3j we translate 
these results into implied volatility asymptotics and prove formulae (JT|) and © above. In Section 2] we 
calibrate the Heston model and provide numerical examples based on formulae ([TJ and ([2]). 

2 Large-time behaviour of call options 

Throughout this article, we work on a model (0, J-, P) with a filtration (J r t)t>o supporting two Brownian 
motions, and satisfying the usual conditions. Let (St)t>o denote a stock price process and we let X t := 
log(Sf). Interest rates and dividends are considered null. We assume the following Heston dynamics for 
the log-stock price: 

dX t = -\Y t dt + VY t dW t \ X = .to e R, 

dY t = k (6 - Y t ) At + aVYdW?, Y = y > 0, (3) 
d(W\W 2 ) t = pdt 

with k, 9, cr, i/o > 0, \p\ < 1. 

The Feller condition 2k9 > a 2 ensures that is an unattainable boundary for the process Y. If this 
condition is violated zero is an attainable, regular and reflecting boundary (see chapter 15 in [22] for the 
classification of boundary points of one-dimensional diffusions). Since the analysis in this paper relies 
solely on the study of the behaviour of the Laplace transform of the process X, which remains well defined 
even if the Feller condition is not satisfied, we do not assume that the inequality 2k8 > a 2 holds. 

Let us now define S := k — pa, p := yl — p 2 , and 8 := kO/R. Throughout the whole paper, we will 
assume R > 0. This assumption ensures (see Theorem 2.1 in |12j ) that moments of S greater than 1 exist 
for all times t. This condition is fundamental for the analysis in the paper and is usually assumed in the 
literature (see [53] and [3J). When this condition is violated the limiting logarithmic Laplace transform 
V defined in (j4} of the process S does not have the same properties, and further research is needed to 
understand how the implied volatility behaves in this case. We know from [3J and [12] that R is the 
mean-reversion level of the process Y under the so-called Share measure that is equivalent to the original 
probability measure with the Radon-Nikodym derivative given by the share price itself. If R < 0, the 
process Y will be neither mean-reverting nor ergodic under the Share measure. In the equities market 
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this does not constitute a problem since the calibrated correlation p is always negative. However this 
assumption may be restrictive in markets such as foreign exchange, and further research is required to 
relax it. Let V be the limiting logarithmic moment generating function of X defined as 

V(p) := lim r 1 logE (cxp (p (X t - x ))) , (4) 



for all p such that the limit exists and is finite. It follows from [5] that V is a well defined and strictly 
convex function on (p_,p + ) and is infinite outside, where 



p± := i-2np + a± v/cr 2 + An 2 - Anpaj / (2a -f) , (5) 
with p_ < and p+ > 1. Furthermore the function V takes the following form 

V(p) = ^ (k - pap - d(-ip)j , iov p G (p-,p+) , (6) 

where 

d(k) := y/( K - ipak) 2 + a 2 {ik + k 2 ), for k G C, (7) 

and we take the principal branch for the complex square root function in ([?]) ■ 

Let us now define the Fenchel-Legendre transform V*(x) := supjpx — V(p) : p G of V, 

which was computed in |12j and is given by the formula 

V*(x) =p*(x)x-V(p*(x)), forallxGR, (8) 

where the function p* : R — > (p_ ,p+) is defined by 

1/2 

a - 2np+ (n6p + xa) ({a 2 + An 2 - Anpa)/(x 2 a 2 + 2xn6pa + k 2 9 2 )) ' 

~2a~p 2 

Tedious but straightforward calculations using the explicit formulae above yield the following proposition. 



p*(x):= — - — o __ 2 " x ^ , forallxGR. (9) 



Proposition 2.1. The function p* : R — > (p—,p+), where p- and p+ are defined in l[5]). is strictly 
increasing, infinitely differentiable and satisfies the following properties 

p*(-9/2)=0, p*(e/2)=l, lim p*{x)=p- and lim p*{x)=p + , 

x — V — oo x— f+oo 

as well as the equation 

V'(p*(x))=x, forallxeR. (10) 

Since the image of p* is (p_ , p + ) , the function is well defined on R. The following properties of 
V* arc easy to prove and will be used throughout the paper: 

(a) V*'(x) = p*(x) for all x G R; 

(b) F*" (x) > for all x G R; 

(c) x n- is non-negative, has a unique minimum at — 9/2 and V*{— 9/2) = 0; 
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(d) x i— > V*(x) — x is non-negative, has a unique minimum at 0/2 and V* (0/2) = 0/2. 

From the definition ([5]) of V™ and relation (fTU|) . the equality in (a) follows. The inequality in (b) is a 
consequence of (a) and Proposition 12.11 Now, (a), (b) and Proposition 12.11 imply that — 0/2 is the only 
local minimum of the function V* and is therefore a global minimum. The definition of V* given in ([5]) 
implies V* (—6/2) = —V(0) = 0. Since the stock price S is a true martingale (see [3]), we have V(l) = 
and Proposition 12.11 implies that V* (0/2) — 6/2 > 0. This proves (c). From (a) and Proposition we 
know that the function x i-> V*(x) — x has a unique minimum attained at 6/2 and V* (6/2) —6/2 = 0. 
Therefore (b) implies (d). 

2.1 Large-time behaviour of call options under the Heston model 

In this section, we derive the asymptotic behaviour of call option prices under the Heston dynamics ([3]) 
as the maturity t tends to infinity, both for maturity-dependent and for fixed strikes. Before diving into 
the core of this paper, let us introduce the function I:Kx R + x R 2 -> 1 by 

1 (x, t; a, b) := (l - c xt ) l {x<a} + l {a<x<b } + \\x=b] + (l - e a *) l {x=a }, (11) 



which will feature in the main formulae of Theorem 12.21 and Proposition 12.71 The next theorem is the 
main result of the paper and its proof is given in Section [5j 

Theorem 2.2. For the Heston model and the assumptions above, we have the following asymptotic 
behaviour for the price of a call option with strike Sq exp(xi) for all 



■J-E (St - S e xt ) + = X (x, t; -6/2, 6/ 2) +(2^)^' 2 exp i-(V*(x) - x) t)A(x) (l + O (l/t)) , as t 



00, 



where 

U(p*(x)) 



A(x) 



y/V"(jf(x)) 



if xeR\ {-6/2,6/2} , 



p*(x) (p*(x)-iy 
-l-sgn(x) QQ^-tfVW)) , ifxe {-8/2,0/2} 



(12) 



where 

V is defined in p* in @, V* in ((SJ) and d in ([7]) and sgn(x) equals 1 if x is positive and —1 otherwise. 



Remark 2.3. Property (b) on page 2] implies that the square root \/V" (p*(x)) is a strictly positive real 
number. Note also that the function A defined in (fT2|) is not continuous at the points —0/2 and 0/2 by 
Proposition ^. II (see also Figured] on page ITTj) . 

Remark 2.4. It is proved in Proposition [231 that p* (-0/2) = and p* (6/2) = 1. Note further that 
U'(0) = (6- y )/(2«) and U'(l) = (y - 6)/(2h). 
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Remark 2.5. Theorem l2.2l is similar in spirit to the saddlepoint approximation for a density of random 
variable X given in Butler [8] 

f x (x) « (2nK"{p*{x))) ^ cxp (K(p*(x)) - xp*{xj) , 

where K(p) := logE(exp(pX)), and p*(x) is the unique solution to K'(p*(x)) = x. Here X := X t — xo, 
-K"(p) = tF(p) + 0(1), and we substitute x to xt so that 

fc- I0 (xt) a (2 7 rF"(p*(x))t)- 1/2 cxp(F(p*(.T))t-xp*( a; )t) - (27TF" (p*(x))i)" 1/2 cxp(-^* (x)t). 

In order to precisely compare our result to the existing literature, we prove the following lemma, which 
gives the asymptotic behaviour of vanilla call options when the strike K is fixed, independent of the 
maturity. The following lemma was derived in [27], Chapter 6; a rigorous proof is detailed here in 
Appendix DO 

Lemma 2.6. Under the same assumptions as Theorem \2.2l for any x G R, we have the following 
behaviour for a call option with fixed strike K = Sq exp(x), 

■^-E (St — K) + = l + 42!eacp((l-p*(0))z-V*(0)t) (l + 0(l/t)), as t -> oo. 

2.2 Large-time behaviour of the the Black-Scholes call option formula 

By a similar analysis, we can deduce the large-time asymptotic call price for the Black-Scholes model. 
This result is of fundamental importance for us since it will allow us to compute the implied volatility 
by comparing the Black-Scholes and the Heston call option prices. Throughout the rest of the paper, we 
let Cbs(<So> K> tr denote the Black-Scholes price of a European call option written on a reference stock 
price S, with strike K > 0, initial stock price So > 0, timc-to-maturity t > 0, and volatility £ > (with 
zero interest rate and zero dividend). Similar to Section [2] let us define the function Vbs : R — !> R as in 
(J3J, where X t := log(St). In the Black-Scholes case, it reads 

Vbs(p) =P(P- l)E 2 /2, for all pGM. (14) 

Similarly to flSJ, we can define the functions V^ s : R x R^ — > R and j»g S : R — > R, by 

V£ s (x, S) := (x + S 2 /2) 2 / (2S 2 ) , for all £ G R% (15) 

and 

PbsO) : = (x + £ 2 /2) /£ 2 , for all x G R. (16) 
The following proposition, proved in Appendix [Bj gives the behaviour of the Black-Scholes price as the 
maturity tends to infinity. 

Proposition 2.7. Let a > and let the real number a\ satisfy a\ > —a 2 t for large times t. Then for 
all i £ 1, toe have the following asymptotic behaviour for the Black-Scholes call option formula in the 
large-strike, large-time case 

±-C BS (s ,Soe*V,y^) = X (x, t; - g , g ) + ^ e~ W^)->«(1 + 0(1/*)), 



G 



where 



A B s{x,a,ai) := exp ( -on ( — r - 1 ) ) j^l^/i^/a} h ; l{ x =±a*/2}, (17) 



and where the function I is defined in Qll[) . 

Remark 2.8. If we set a\ = in Proposition 12.71 we obtain the large-time expansion for a call option 
under the standard Black-Scholcs model with volatility a and log-moneyness equal to xt. 

As in the Heston model above, we derive here the equivalent of Proposition 12 . 71 when the strike does not 
depend on the maturity anymore. 

Lemma 2.9. With the assumptions above, we have the following behaviour for the Black-Scholes call 
option formula in the fixed-strike, large-time case 

-^C BS (Sq, S c x , t, ^<7 2 + ai /t) = 1 - -^2= exp (-<r 2 t/8 + x/2 - oi/8) (1 + 0(1/ 1)) . 
dq \ / ay/lTt 

This lemma is immediate from the Black-Scholes formula given in Appendix [B] and the approximation 
(|A-2[) for the Gaussian cumulative distribution function. 

3 Large-time behaviour of implied volatility 

The previous section dealt with large-time asymptotics for call option prices. In this section, we translate 
these results into asymptotics for the implied volatility. Recall that [12] and [27] have already derived 
the leading order term for the implied volatility in the large-time, fixed-strike case. Our goal here is 
to obtain the leading order and the correction term in the large-time, large-strike case. Theorem 13.11 
provides the main result, i.e. the large-time behaviour of the implied volatility in the large strike case. In 
the following, a t (x) will denote the implied volatility corresponding to a vanilla call option with maturity 
t and (maturity-dependent) strike 5*0 exp(xi) in the Heston model ([3]). We now define the functions 
&%, : R ->• M + and Si : K -> R by 

*Ux) := 2 (2V*(x) -x + 2 (l xe( _ e/2m - l xm \(-e/2,e/2)) VV*(x) 2 - V*(x)x) , for all xel, (18) 

and 

(x'/atix) - 1/4J log {A(x)/A BS (x, &oo(x),0)), x€R\ 



2' 2 



ai(x) := 2 < 



Zoo(x) / . , (V"'{p*{x)) 



(19) 



y/V"{p*{x)) V \6V"(p*(x)) \ 2 ' 2. 

where A is defined in (jUJ, A BS in (pi]), U in $T5fy, V in ([6]), p* in ^ and in ©. They are all 
completely explicit, so that the functions ct^ and Si are also explicit. From the properties of V* proved 
on page|H V*(x) and V*(x) — x arc non- negative, so that a 2 Ki (x) is a well defined real number for all 
x £ R. Then the following theorem holds. 



7 



Theorem 3.1. The functions and a\ are continuous on K and 

a t [x) = (x) + a± (x) jt + o (1/t) , for all x G R, as t oo. 

Furthermore the error term \af (x) — ^^(x) — a\(x)/t\t tends to zero as t goes to infinity uniformly in x 
on compact subsets ofM\{ — 9/2, 9/2}. 

Proof. We first prove that the functions er^ and a\ are continuous. In fact, the continuity of the function 
(Too follows from properties (c) and (d) on page 2] We have already observed that the function A defined 
in (|12l) is discontinuous at 9/2 and —9/2 (see Remark I2.3[) . Elementary calculations show that in the 
neighbourhood of 9/2 we have the following expansion (where := (9 /V" (l)) 1 ^ 2 ) 

al(x) = + 2 (1 - 9) (x - 9/2) + ^- (l - I + ^^e) (x - 9/2) 2 +0(\x- 9/2\ 3 ) . 

This expansion, together with formula (TTT]) . implies that the function x h-> Abs (x, CToo(^), 0) is also 
discontinuous at 9/2 and —9/2. Since the quotient x A(x)/A-Bs{x,6^ aD (x),0) is strictly positive on 
R \ {—0/2,0/2}, the function ai is therefore continuous on this set. The expansion above can also be 
used to obtain the following expression 

A , A{X \ n^ = 1+ T7^(^ 1 )- 1 -SrT + ^ N ) (--^72) + O(|x-0/2| 2 ), 

-4Bs(^,croo(a;),0) V"(l) \ 6V"{1) 6/ 

which implies the equality lim ai(x) = ci\ (9/2). A similar argument shows continuity of Si at —9/2. 

z-s-e/2 

We now prove the formula in the theorem in the case x > 9/2 (recall that 9 > 0). Note that (Too (a:) 
as defined in (|18|) satisfies the following quadratic equation 

V*(x)-x = V£ s (x,a 00 (x))-x, for all x el, (20) 

where V* is given by ((HJ and V^ s by (fT5|) . The proof of the theorem is in two steps: first, we prove 
the uniform convergence (on small intervals) of the implied variance to a 2 K> (x), then we derive a similar 
result for the first order correction term. As a first step, we have to prove that, for all S > 0, there 
exists t*(5) > such that for all t > t*(5) and all x in some small neighbourhood of x 1 > 9/2 we have 
|<Tt(x) — (Too (a;) | < 8. By Theorem 12.21 and (|20|) we know that for all e > 0, there exists t*(e) and an 
interval containing x' such that for all x in this interval and all t > t*(e) we have the lower bound 

cxp ( - (V£ s (x, a^x)) - x + e) t) = exp ( - (V*(x) - x + e) t) < i-E (S t - S e xt ) + , (21) 
and the upper bound 

jrE (S t - S e xt ) + < exp ( - (V*{x) - x - e) t) = exp ( - {V£ s {x, a^x)) -x-e)t). (22) 
Note that 

ct^(x) - 2x = 4(V*(x) - z) - VO^fc) - x) 2 + (V*(x) - x)x} < 0, for all z G M\(-6/2, 6/2), 
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since V*(x) — x > by property (d) on page 2] For x fixed, the function S H> V£ s (x, S) — x denned on 
(0, V2x) is continuous and strictly decreasing, where V£ s is given in (|15p . Thus, for any S > such that 
<7oo (a;) ± i5 £ (0, y/2x) for all x in the small neighbourhood of x' define 

e 1 (8) := (V B * S (x, (ar)) - V£ s (x, (x) + 5))/2 > 0, 

£2(<5) := (v£ s (x,a oa (x)-5)-V£ s (x,& oo (x)))/2>0. 

Note that the limits lime! (<$) = lim e 2 (<$) = are uniform in x on the chosen neighbourhood of x' . 

8— s-0 S— >0 

Combining (j2"Tj) . (|2"21 and Proposition 12.71 there exists t*(S) such that for all t > t*(6) and all x near a/ 
we have 

t^Cbs (5 , Sbe", t, avoix) - 5) < exp ( - (V^z, a*, (a) - S) - x - e 2 (S))t) < -^E (S t - S e xt ) + , 
and 

-^E (S t - S c xt ) + < exp ( - (Vggfoo-ooCs) +5) - a: + ei(a))f) < -^C BS (Sb.Soe**, ^M*) + 6) . 

DO V /Do 

Thus, by the monotonicity of the Black-Scholes call option formula as a function of the volatility, we have 
the following bounds for the implied volatility a t {x) at maturity t 

0-00(2;) - 5 < a t {x) < &oo(x) + 5. 

In the second step of the proof we show that for all 6 > there exist t* (5) and a small neighbourhood 
of x' such that, for all t > t*(S) and all x in this neighbourhood, the following holds 

\a 2 (x)-a 2 OQ (x)-a 1 (x)/t\<S/t. (23) 

Note that the definition in (|T9"| implies 

A(x) = ABs(x,(Toc(x),di(x)), for all 

Theorem 12.21 implies that for any e > there exist a small interval containing x' and t* (e) such that for 
all x in this interval and all t > i*(e) we have 

— E (S t - S Q c xt ) + < Abs{x ' e -(yBs(*.*~(*))-*) y, (24) 

So \j2nt 

and 

fjgs^ 1 £oo(g) 1 ai^)) e -(vg 8 (x,g ae (x))- 8 )t e - e < ^ E ( 5t _ 5- oe **) + 

Let 5 > and define e(S) := (4x 2 /&^ (x) — l) (5/16 > 0. The reason for this definition lies in the following 
identity 

A B s(x,o- 00 (x),a 1 (x))e ±e(,5) = A BS (x, a^x), a^x) ± S)c Te{S) for all x6l 

which holds by definition ([T7|) . By (|24|) there exist i * (<5) and a neighbourhood of x' such that for all x in 
this neighbourhood and all t > t* (<5) we have 

— E fSt - 5 , c a;t ) + < ^bs(z, 0-00(2:), ai(x) + j) e -(y B ' 3 ( a; ,fro ( a ;))- a! )t e -6(g) 
So V y V^rt 



< ^C BS (s , S e xt , t, y/aUx) + (ai(x)+S)/tj 
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and 

S V ' ^/2^rt 

> ^<?bs (s ,S e xt ,t, ^/&Ux) + (ai(x)-S)/t) . 

By strict monotonicity of the Black-Scholcs price as a function of the volatility, we obtain 
<^o(z) + (&i(s) - S)/t < a 2 (x) < 8*,(x) + (ai(x) + 8)/t, 

for all x in some interval containing x' and all t > t*{S). This proves ([23]). An analogous argument 
implies the inequality in (|23p for x' <G M\{—6/2, 0/2} since the functions A and a; i— > Abs(x, a 00 {x),ai(x)) 
are continuous on this set and by ([2"U|) the following identity holds 

I(x,t;-6/2,6/2)=I(x,t;-&l (x)/2,al (x)/2) for all (i,f)elxl + 

(the function X is defined in (fTTjO . In the cases x' £ {—0/2,0/2} a similar argument proves (f2"3"]l at the 
point x' but not necessarily on its neighbourhood. 

All that is left to prove is that the inequality in ([23]) holds uniformly on compact subsets of the 
complement R\{ — 0/2, 9/2}. To every point in a compact set we can associate a small interval that 
contains it such that inequality ([23]) holds for all x in that interval and all large times t. This defines a 
cover of the compact set. We can therefore find a finite collection of such intervals that also covers our 
compact set. It now follows that the inequality in ((23)) holds for any x in the compact set and all times 
t that are larger than the maximum of the finite number of t* (8) that correspond to the intervals in the 
finite family that covers the original set. □ 



3.1 The large-time, fixed-strike case 

This section is the translation of Lemmas 12.61 and 12.91 in terms of implied volatility asymptotics, and 
improves the understanding of the behaviour of the Heston implied volatility in the long term. Let 
<Jt{x) denote the implied volatility corresponding to a vanilla call option with maturity t and fixed strike 
K = So exp(x) in the Heston model ([3]). Let us define the function a\ : M. — > R by 

diO) := -8 log (-A(0)y/2V*(of) +4(2p*(0) - l)x, for all x G R, (25) 

where A is defined in (IT21) . V* in ([8]) and p* in . Elementary calculations show that A(0) < 0. From 
the properties of V* on pagelH a\(x) is then well defined as a real number for all i£R. 

Theorem 3.2. With the assumptions above, we have the following behaviour for the implied volatility in 
the fixed-strike case 

a\ (x) = 8V*(0) +a 1 (x)/t + o{l/t) , for all x G R, as t -> oo, 

where V* is given by ([SJ and a\ by (|25[) . The error term \a 2 t {x) - 8F*(0) - a\(x)/t\t tends to zero as t 
goes to infinity uniformly on compact subsets of R. 
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Figure 1: We plot the two functions A (left) and Si (right) with the calibrated parameters on page IT21 
The two crosses on the left graph represent the values A(— 0/2) and A(0/2) at —8/2 w —0.025 and 
0/2 ps 0.022 where the function A is discontinuous. 



Proof. The proof of this theorem is similar to the proof of Theorem 13.11 It is in fact simpler because we 
do not have to consider special cases for x' € M. We therefore only give an argument that implies the 
inequality 

cr t 2 (x) - 8V* (0) - oi (a;) /t<6/t (26) 

holds on some small neighbourhood of x' for all large t. Lemma T2.6I and (|2"S"j) imply that for all e > 0, 
there exists a i*(e) and a neighbourhood of x' such that for all x in this neighbourhood and all t > t*(e) 
we have 

(5 t - S* e x ) + < f 1 + (2^t)- 1 / 2 A(0) exp (x(l - p*(0)) - V*(0)t)) e e 
= | 1 1 exprx^-^f^/S-^rOli) I e e . 

\ 2^7T^*(0)i 7 

For any 8 > a continuity argument implies the existence of e(<5) > 0, which tends to zero uniformly on 
a neighbourhood of x', such that 



^/2-a 1 (x)/8-y*(0)t „e(*) 



2^7ry*(0)i 



=e x/2-( ai (x)+5)/8-V*(0)t e -e(5)_ 



Therefore there exists t*(<5) > and an interval containing x' such that for all x in this interval and all 
t > t*(S) we have 

(St - c x ) + < -trC BS (So, S Q c x , t, y/8V*(0) + ( ai (x) + 6)/t) . 

The monotonicity of the Black-Scholcs formula as a function of the volatility implies (f2"6")l . The proof can 
now be concluded in the same way as in Theorem 13.11 □ 
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4 Numerical results 



We present here some numerical evidence of the validity of the asymptotic formula for the implied 
volatility obtained in Theorem 13.11 We calibrated the Heston model on the European vanilla options 
on the Eurostoxx 50 on February, 15th, 2006. The maturities range from one year to nine years, the 
strikes from 1460 up to 7300, and the initial spot So equals 3729.79. The calibration, performed using the 
Zcliadc Quant Framework, by Zeliadc Systems, on the whole implied volatility surface gives the following 
parameters: « = 1.7609, 6 = 0.0494, a = 0.4086, y = 0.0464 and p = -0.5195. We then use these 
calibrated parameters to generate an implied volatility smile for two different maturities, 5 and 9 years 
(note that for convenience the smiles were generated with interest rates and dividends equal to zero). 
The two plots below contain the generated implied volatility smile for the Heston model for the two 
maturities as well as the zeroth and first order approximations obtained in (|18p and in Theorem 13.11 
The errors plotted are the differences between the generated Heston implied volatility and the implied 
volatility obtained by the two asymptotic approximations. 



5 Proof of Theorem [23] 

The proof of the theorem is divided into a series of steps: we first write the Heston call price in terms 
of an inverse Fourier transform of the characteristic function of the stock price ([27]) . Then we prove a 
large-time estimate for the characteristic function (Lemma (5Tj. The next step is to deform the contour 
of integration of the inverse Fourier transform through the saddlepoint of the integrand (Equation ([9]) and 
Proposition ^. ip . Finally, studying the behaviour of the integral around this saddlepoint (Proposition l5.8[) 
and bounding the remaining terms (Lemma I5.6[) completes the proof. The special cases x = —8/2 and 
x = 8/2 in formula (fT2j) are proved in Sections 15.31 and 15.41 

5.1 The Lee-Fourier inversion formula for call options 

Using similar notation to Lee [26], set 

A t:X ■= [v G K : E (exp (v(X t - z ))) < oo j , for all t > 0, 

and define the characteristic function <j>t : M — > C of Xt — xo by 

<f>t{z) := E (exp (iz(X t - a- ))) , for all t > 0. 

From Theorem 1 in |28j and Proposition 3.1 in [3J, we know that (f>t(z) can be analytically extended 
for any z € C such that — 3(z) € A tt x, and from our assumptions on the parameters in Section [51 
(p_,p_l_) C At t x for all t > 0. By Theorem 5.1 in [26], for any a £ (p-,p+), we have the following Fourier 
inversion formula for the price of a call option on St 

^E (S t - K) + = <M-i)l { o<«<i} + 0M- 1 ) - c"0t(O))l { « <o} + -0 4 (-i)l {a=1} 
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Figure 2: The left plot represents the leading order term (Too (dashed) defined in ([18]). the asymptotic 
formula in Theorem 13.11 (solid) and the true implied volatility (crosses) as functions of the strike K for 
a maturity equal to 5 years. On the right, we plot the corresponding errors between the true implied 
volatility and (dashed) and between the true implied volatility and our formula (solid) . The parameter 
values are given Section S) 
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Figure 3: The left plot represents the leading order term &oo (dashed) defined in (fl8|) . the asymptotic 
formula in Theorem 13.11 (solid) and the true implied volatility (crosses) as functions of the strike K for 
a maturity equal to 9 years. On the right, we plot the corresponding errors between the true implied 
volatility and (dashed) and between the true implied volatility and our formula (solid). The parameter 
values are given Section |U 



13 



where x := \og(K / So) and 7+ : M — > C is a contour such that 7+(w) := u — i(a — 1). The first four terms 
on the right hand side arc complex residues that arise when we cross the poles of (iz — z 2 ) at z = 
and z = i. We now set k = i — z, substitute x to xt, and use the fact that St is a true martingale for all 
t > (see Proposition 2.5 in [3]). From now on, as k will always denote a complex number, we use the 
notation k = k r + iki for k r , fcj G R. Note that 



±k — k 2 J \ \ ±k — k 2 

that k r H> ^ c lfc:rt p J is an even function and k r M> 3 ^ c ifcz* e -i-l^ Q> ^ an dd function. 

Clearly the normalised call price S , _1 E (St — So exp(xt)) + is real, so if we take the real part of both sides 
and break up the integral, we obtain 

j-E (S t - Soe^f = 1 {0<q <i } + (1 - e"*) 1 {q<0} + ±l {a=1} + ±eA l {a=0} 
, exp{xt) ( ( f f \, kxt M- k ) 



2» Hlt + /J e S^ d *J- (27) 

for any R > 0, where, for any a £ K, we define the contours 

7c : (—00, — R] U [R, +00) — > C such that 7 a (u) := u + ia, (28) 

and 

Ca : (-R, R)^C such that Ca(«) := u + ia. (29) 

For ease of notation, we do not write explicitly the dependence of these contours on R. We will see later 
how to choose R. In the following lemma, we characterise the large-time asymptotic behaviour of the 
characteristic function <p t . 

Lemma 5.1. For all k £ C such that —hi £ (p_,p + ) ; we have 

4>t(k) = cxp (v(ik)tj U(ik) (1 + e(fc, t)) , as t -> 00, 

ft{d{k)) > 0, and e{k,t) = O ( e - tSR ( d ( fe ))), wh ere U is defined in (fT3"|). p + in (|5|), d in ([7]) and V is 
the analytic continuation of formula ([6]). 

If — fcj is not in (p_,p + ), this large time behaviour of <pt still holds, but ^R(d(k)) might be null (for instance 
if k r = 0) so that e(k, t) does not tend to as t —> 00. 

Proof. From [2], we have, for all k £ C such that —hi £ (p_,p + ), 

«*> — - (^T)) - (s-^'i 1 ^) • <*» 

where <i is defined in ([7]), V is the analytic extension of formula (O and the correct branch for the 
complex logarithm and the complex square root function is the principal branch (see also [2] and |28j ) 
and g : C — > C is defined by 

g( k ) := K ~ lpal !~tnh forallfceC. (31) 
w K- ±pak + d{k) y ' 
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For all k G C such that — fcj G we have 3?(d(fc)) > 0. Let 



(k — pcrifc + 
Then wc have 

f (fc) =exp ^(ifc)^t/(i/c)ei(fc,<)e 2 (fc,<): for all £ > 0. 
Then, as t tends to infinity we have 

C1 (M) = 1 + 9 =^e~ d ^ + O (c- 2 *)') and e 2 (k,t) = 1 + _ i + O f(e d «< - l)~ 2 ) 

for some constant c. Set e(k,t) := ei(k,t)e2(k,t) — 1 and the lemma follows. □ 

5.2 The saddlepoint and its properties 

Wc first recall the definition of a saddlepoint in the complex plane (see [7]): 

Definition 5.2. Let F : Z — > C be an analytic complex function on an open set Z. A point zq G Z such 

dF 
dz 



that the complex derivative 4^ vanishes is called a saddlepoint 



Note that the function V : — > M defined in ((4]) can be analytically extended and we define for 
every id the function F x : Z — > C by 

F x (k) :=-±kx-V(-±k), where Z := {k £ C : h £ {p-,p+)} ■ (32) 

Note that the exponent of the integrand in (|2T|) has the form —F x (k)t by Lemma 15.11 therefore the 
saddlepoint properties of F x given in the following elementary lemma are fundamental. 

Lemma 5.3. The saddlepoints of the complex function F x : Z — > C are given by 
± a - 2np± (n9p + xa)n(x 2 a 2 + 2xK0pa + k 2 ^ 2 )" 1 / 2 



where r\ := \J a 2 + 4k 2 — Anpa. 

Proof. Since we are looking for the saddlepoint in Z, we can use the representation in © for the function 
V. Therefore the equation F' x (z) — is quadratic and hence has the two purely imaginary solutions 
Zq(x) since the expression (x 2 a 2 + 2xnQpa + k 2 6 2 ) = (xa + nOp) 2 + n 2 8 2 (l — p 2 ) is strictly positive for 
any x G R. It is also clear from the definition of p- and p+ given in ([5]) and the assumptions on the 
coefficients that 3 (^(x)) G and therefore Zq (x) G Z are saddlepoints of F x . □ 

The next task is to choose the saddlepoint of the function F x in such a way that it converges to the 
saddlepoint of the function F x s (k) := —ikx — Vbs(— ifc) for all k G Z, where Vbs is given by (fl"4"]). in 
the Black-Scholcs model as both the volatility of volatility and the correlation in model (|3|) tend to zero. 
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It is easy to see that the saddlepoint of F® s equals ipQ S (x) for any were p^ s is given by (|TB)) . Wc 

can rewrite 3 (z (a;)) defined in Lemma 15.31 as 

^, +. 1 ko , K0pn(x 2 a 2 + 2xK9pa + ^e 2 )- 1 / 2 , xan(x 2 a 2 + 2xn6pa + k 2 B 2 Y x I 2 
^(^) = V-^ ± " 2aJ 2 ^ >—, (33) 

where p is defined page [3] The first term converges to 1/2 and the last one to ±x/9 as (p,a) tends to 0. 
When both p and a tend to 0, a Taylor expansion at first order of the third term gives 

n6pr](x 2 a 2 + 2xn9pa + ^B 2 )- 1 ' 2 pr\ 



2a p 2 2a p 2 ' 



(34) 



Take now the positive sign in p3[) . then the second and third terms cancel out in the limit because r\ 
converges to 2k as (p, a) tends to 0. In that case we have lim 9 (zq(x)) — 1/2 + x/6 — Pbs( x ) f° r a ^ 

(p,<r)^0 

x G K, where the Black-Scholes variance is equal to 9. If we were to take 3 (zq(x)) for the saddlepoint. 
in the limit we would not recover p^g(x) by (|34|) . since the function (p,a) > p/a has no limit as the pair 
(p, a) tends to 0. Therefore we define the saddlepoint to be Zq(x). Moreover we observe the following 
equality 

$S(z£(x))=p*(x), for all x el, 

where p*(x) is defined in ©. 

Remark 5.4. In the fixed strike case, i.e. when x = 0, we obtain 

-2np + a + p7j 
P (0) = 2aJ 2 ' 

The corresponding saddlepoint ip*(0) is the same as the one in Chapter 6 in |27j . 

The following lemma is of fundamental importance and will be the key tool for Proposition 15.81 

Lemma 5.5. Let k G Z. Then, for any fcj G (p_,p + ), the function k r \— ¥ 5i (— ikx — V(— i-k)) has a 
unique minimum at and is strictly decreasing (resp. increasing) for k r G (— oo,0) (resp. k r G (0,oo)J. 

Proof. Note that the statement in the lemma is equivalent to the map k r i— > — (V(— i(A; r + ifcj))) having 
a unique minimum at k r = for any fej G and being increasing (resp. decreasing) on the positive 

(resp. negative) halfline. Let fcj G then 

SR(V(-i(fc r + ifc 4 ))) = 4 ( K - ^ - K (vM^) + i«(M)) , 

where 

u{k r ) := a 2 p 2 k 2 — a 2 p 2 k 2 — a(2np — a)ki + k 2 and v(k r ) := (2np — a + 2ap 2 k^ ak r . 
From the identity and the fact that the principal value of the square-root is used, we get 

3? (yju{k r ) + iu(£v)) = iy / 2u(fc r ) + 2y/u 2 (k r ) + v 2 (k r ) (35) 



16 



is monotonically increasing in u, u 2 and v 2 . First, note that u 1 (k r ) = 2o 2 p 2 k r , hence u is a parabola with 

a unique minimum at k r = 0, so that, from (|35[) . it suffices to prove the following claim: 

Claim: For every k% G (p—,p+), the function g := u 2 + v has a unique (strictly positive) minimum 

attained at k r = and is strictly increasing (resp. decreasing) for k r > (resp. /c r < 0). 

Let us write u(k r ) = o 2 p 2 k 2 + ip(ki), for all k r G R, where ^(fc,) := k 2 — o 2 p 2 k 2 — o(2np — a)ki. We have 

5 (fc r ) = o 4 p 4 k 4 + (2p 2 iP{h) + (2np - a + 2crp 2 fc i ) 2 ) a 2 k 2 r + ^(h) 2 for all k r G R. (36) 

The coefficient cr 2 p 4 and the constant k 2 are strictly positive, so the claim follows if x(k%) > f° r a ^ 
fc, G (p_,p + ), where 

X(fci) := (2p 2 ip(ki) + (2np -a + 2op 2 ki) 2 ^j = 2o 2 p 4 k 2 + 2ap 2 (2np - a) h + k 2 + (2np - of . 

The discriminant is A x = — Ao 2 p 4 {2k 2 + (2np — er) 2 ^ < 0, so that \ nas no rca l r0 °t an d is hence always 
strictly positive. This proves the claim and concludes the proof of the lemma. □ 



The following two results complete the proof of Theorem 12. 2\ by studying the behaviour of the two 
integrals in (|27[) as the time to maturity tends to infinity. The following lemma proves that the integral 
along "fp*( x ) is negligible and Proposition 15.81 hereafter provides the asymptotic behaviour of the integral 
along the contour Cp*(x)- 

Lemma 5.6. For any i£l and any m > V*{x), there exists R(m) > such that for every k G Z with 
\k r \ > R(m), we have 

|exp (ikxt) 4>t{— k)\ < exp(— mi), for allt>\. (37) 



Therefore 



where the contour r y p -( x ) is defined in ([28 



1 1 e ^M^) di 



O f e -( m - x A , (38) 



Remark 5.7. (i) For every x G R, we have V*(x) > x by (d) on page|3]and hence m — x > 0. Therefore 
the modulus of the integral (|38[) tends to zero exponentially in time and in m. 

(ii) Recall that p*(x) G by Proposition 12.11 and hence inequality (|37)) can be applied when 

estimating integral 



Proof. We only need to prove (|3"8|) . Recall from Lemma I5TT1 after some rearrangements, that 
fc(-fc) = e (t+y /(ne))v(- ik ) (1 _ g( _ fc)) ^/- 2 ^ + Q ^-m W - k)) ^ 

It follows from equations (J7J), © and (j3"T|) that 

5ft (d (— (k r + ifcj))) ~ crp ! as —> oo, 
V (— i(fc r + ifcj)) ~ —KOp \k r \ /o, as |fc r | — > oo, 
lim <? (— (k r + iki)) = (p — lp) 2 ^ 1, since \p\ < 1. 
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Hence there exists a constant C > 0, independent of fc and t, such that the following inequality holds 

\c ikxt (f)t(-k)\ < Cexp(-t(kiX- I + a' 1 n6p\k r \)y 

Define R(m) := max{<7 (m + 1 — kiX + log(C)) /(n0p), 1}. Then if |fc r | > max {R(m), R}, where the 
positive constant R is given in definition ([2"5]l . the equality ([55)) follows. □ 



Proposition 5.8. For any R > tme? i6l\{-S/2,9/2}, iue Ziawe as t oo, 

exp(rf) -5R ( / e^#4ld^ = ^(-(V^)-*)*) (A(a;) + o {l/t)) f m 



2tt i ftn ifc -A: 2 / 72?rt 



where A is given in (|12[) . T^* m ([8J and £p*(z) ^ n (|29j) . 

Proof. Let \ {— 0/2,0/2}. Applying Lemma [5~T1 on the compact interval [— we have 

e ibf M dfc = / ^h^ c (ife-+v(-ifc))t (l + e (jfc, t)) dfc, 

for i large enough. By Lemma [5751 we know that fc r i-> — 5R(i(fc r + ip*(x))x + V^(— i(fc r + ip*(x)))) has a 
unique minimum at fc r = and the value of the function at this minimum equals V*(x) by the definition 
of V* . The functions V and {/ are analytic along the contour of integration and thus, by Theorem 7.1, 
section 7, chapter 4 in [53], we have 

/ r u(-±k) (ikx+v( _ ik))t \ £!_ -vw ( u{p*{*)) + 0(1/t) \ 

_«p(-(v(*)-*w (il(a;)+0(1/t)) 



as t tends to infinity. The e(k,t) term is a higher order term which we can ignore at the level we are 
interested in. □ 



Lemma [576] and Proposition 15.81 complete the proof of Theorem 12.21 for the general case. Concerning 
the two special cases, we first introduce a new contour, the path of steepest descent, which represents 
the optimal (in a sense made precise below) path of integration. Note that, the general case can also be 
proved using this path, but Lemma 15.51 simplifies the proof. 



5.3 Construction of the path of steepest descent 

We first recall the definition of the path of steepest descent before computing it explicitly for the Heston 
model in the large time case. 

Definition 5.9. (see [32]) Let z := x + iy, and F : C — >• C be an analytic complex function. 

The steepest descent contour 7 : R — > C is a map such that 

• $t(F) has a minimum at some point zq <E 7 and $t(F"(zo)) > along 7. 
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• $s(F) is constant along 7. 
These two conditions imply that F'(z ) = 0. 

The following lemma computes the path of steepest descent explicitly in the Heston case for the function 
F x given in (|32[) passing through ±p*(x). 

Lemma 5.10. The path of steepest descent 7 in the Heston model is the map 7 : R — > C defined by 

7(5) := s + iki(s), for all s £ R, 

where 



ki{s) := ^ y > ( 4 °) 



/?- (K9p + X(j)y/tp(s) 

13 := k9£(2k P - a), V(s) := 4cr 2 / 5 2 £ 2 s 2 + k 2 9 2 ((2/cp - tr) 2 + 4k 2 ( o 2 ) £ and £ := («0p + xa) 2 + K 2 6» 2 p 2 . 



Note that £ is strictly positive, so that the function fcj is well defined. 

Proof. By definition, the contour of steepest descent is such that the function 3 [F x o 7) remains constant. 
So we look for the map 7 such that Q(F x ('j(s))) = 0, for all s £ R because ^(7(0)) is already real. Using 

/ \-l/2 

the identity 3 (\Ar + iy) = 4 I 2.x + + y 2 ) , for all x 1 y £ R, we find that the function F x o 7 is 

real along the contour 7 : s i— > s + ifcj(s). Note also that this contour is orthogonal to the imaginary axis 
at ki(0) (see Exercise 2, Chapter 8 in [32]). □ 

Remark 5.11. 

• The contour 7 depends on x, but for clarity we do not write this dependence explicitly. 

• The construction of 7 is such that the saddlcpoint defined in ((9]) satisfies ip*{x) = 7(0). 

• We have fcj(s) = 5(7(s)) is an even function of s, i.e. 7 is symmetric around the imaginary axis. 

We now prove Theorem l2.2l in the two special cases x £ {—9/2, 0/2}. In these cases, we need a result 
similar to Proposition l5.8l as Lemma [5751 still holds, i.e. we need the asymptotic behaviour of the integral 
in (|39|) for the two special cases. The problem with these special cases is that (ifc — fc 2 ) -1 in the integrand 
in (|39j) has a pole at the saddlcpoint, so we need to deform the contour using Cauchy's integral theorem 
and take the real part to remove the singularity, before we can use a saddlepoint expansion. 



5.4 Proof of the call price expansion for the special cases 

We here prove Theorem l2.2l in the case x = 9/2 for which p* (0/2) = 1, V* (9/2) = 9/2, and for simplicity 
we also assume that k < (cr — 2p 2 a) /(2p) (the other cases follows similarly). From (|40[) . we see that in 
this case, 7 lies below the horizontal contour 7^ : R — > C such that Jh(s) := s + i (in the other case, 7 lies 
above -fn)- We want to construct a new contour leaving the pole outside. Let e > and 7 e : (— tt, it] —> C 
denote the clockwise oriented circular keyhole contour parameterised as 7 e (0) := i + ec lfl around the 
pole. To leave the pole outside the new contour of integration, we need to follow 7 on R_, switch to the 
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Figure 4: Here we have plotted the closed contour of integration in (|4T|) for e = 0.05, <5 = 0.6, R = 1. 

keyhole contour as soon as we touch it, follow it clockwise (above the pole), and get back to 7 on R + . As 
7 is below 7#, it intersects j e on its lower half, which can be analytically represented as 7" : [— e, e] — > C 
such that 7 e ~( s ) := i + s — i\/e 2 — s 2 . From (|40|). the two contours intersect at s* = ±e. Choose now 
< e < 8 < R ( Lemma 15.121 makes the choice of 8 precise and e must be such that 1 + 2e < p+) and 
define the following contours (they are all considered anticlockwise, sec Figure [4]) 

• Ism. ■ [-R, —5] U [5, R] — > C given by js y R{u) = u + ±ki(S), with ki(S) defined in (|4T)|) ; 

• 7 £! 5 is the restriction of 7 to the union of the intervals [—5, — e] U [e, 8]; 

• 7^ is the portion of the circular keyhole contour 7 e which lies above 7, i.e. the upper half keyhole 
contour as well as the two sections of 7 e between 7 and jh ; 

• ^ are the two vertical strips joining ±i? + i(l + 2e) to ±i? — ±ki(8). 
By Cauchy's integral theorem, we now have 

(L + L + L + L,-L J ^-('w=o, <«, 

Recall that the curves Ci+2 C and 7i+2e are defined in (f2U| and (|2^|) respectively and rewrite (|4*T|) as 




The integral on the left-hand side is equal to the normalised call price 2^^ 1 e et / 2 E \ St — SoC et ^ 2 J by 
Theorem 5.1 in Lee [26], which is independent of e (this holds because 1 + 2e < p + ). For k close to i, we 
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have 



M k) exp(ik9t/2) - |^T+0(i; 



e -«/2 



ifc — fc 2 
so that 

exp (iftflt/2) dfc = (tt + O (e)) c"^/ 2 . (43) 

Lemma f5 . 1 2 1 gives the behaviour of the last integral on the rhs of (|4"2"j) as e tends to for 5 small enough. 

The other integrals can be bounded as follows. By Lemma [5^61 the integral along 7i+2e is O ^e - ^*/ 2 ^ , 
for t > t*(m), R > R(m), as e tends to 0. 

The curves s are both vertical strips of length S and therefore their images are compact sets. 
Applying the tail estimate of Lemma 15.61 along s , we know that for any m > 9/2, there exist t(m) 
and R(m) such that 

/ ^~^l e ikSt/2 dk = O (e- mt ) , for all t > t(m),\k\ > R(m). 
7 r ± ik — k z 

Lemma [5.51 implies that the real function k r n- (— i(k r + iki)9/2 — V(— i(fc r + ifci))) attains its 
global minimum at for any fixed fcj € (p_,p+) and is strictly decreasing (resp. increasing) for k r < 
(resp. k r > 0). It therefore follows that the function u n- sft (—175^(1^)^/2 — V {—i.^s,R{u))\ , where 
u E [-R, —5] U [6, R], attains its minimum value g(8) := 3? ((ki(6) - ±5)9/2 - V{h{5) - iS)), where h(S) 
is defined in (|4"0"|). at the points u = ±<5. It can be checked directly that g(0) = 9/2, g'(0) = and 
g"(0) > and hence for every S > there exists eo > such that the following inequality holds 

9? {-±-is >R {u)9/2 - V(-±y s , R (u))) > 9/2 + e for all u G [-R, -5} U [S, R}. 

Therefore Lemma 1 5. II yields the following inequality 



< e 



-(8/2+e )t 



U(-±k)(l + e(k,t)) 



Ik - k 2 



dk = O (exp (- (9/2 + e ) i)) 



'75, H 

We now prove the following lemma about the integral along j e> s as e tends to 0. 
Lemma 5.12. For 6 > and sufficiently small we have 



^i„t^^ {ikSt/2)ik " \fe"" 2 (- - SWT + H 11 + 0,1/ ' ,) ' 

where U is given by (|13[) and V 6y (|6|). 

Proof. Recall that 7 is the contour of steepest descent defined in Lemma 15.101 and that the curve 7^,5 
is its restriction to the intervals [-5,-e] U [e,8]. Note that s ^ jft ( M^^iM e ij(s)§t/2^ is an even 

function and ( ^iiz2MhL^l e 1 7(s)9t/ 2 ] j s an oc jc{ function. We therefore obtain 

I CX P (^ 2 ) dfc - / * ( ^^f^Tff «P (i7(-)«t/2)) d.s. (44) 

J lc5 lk-k 2 J[-8,-e]u[e,S] V 1 7(s)-7(s) 2 V 
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From (j40|) . for s around 0, we have (i7(s) — 7(s) 2 ) 1 = i/s - 1 + O(s), so 3? ^(i7(s) — 7(s) 2 ) J = 
— 1 + O(s), i.e. taking the real part removes the singularity at k = i. Using Lemma 15.11 we then have 



-<5 



( ^-^ S ^'[% ^Wt/2\ ds = f sjj ( e i 7 ( 8 )e i / 2+ y(-i T(s ))t ds + Q / _ TOt N 

V i7(s) -7(s) 2 J J-s 



for some m > large enough, where we define the function q : M\{0} — > C by 

gw := ^^g))yw forallseM . 

Then, from (|4"0")) . we have the following expansion 

<K S ) = f - - ( + i^C 1 )*) + O (* 3 ) • (45) 



s V6V"(1) 

We can therefore extend the function q to the map q : Bs(0)\{0} — > C for some S > 0, where Bg(0) := {z S 
C : |z| < <5} is an open disc of radius 5. Note that for s E K we have 3? (g(s)) = -1 - V"'{1)/ (6V"(1)) + 
U'(l) + O(s) and hence the function $l(q) : [—6, 6] — > M. does not have a singularity at s = 0. 
Recall that if a function G : Ba(0)\{0} — > C has a Laurent series expansion 



G(z) = ±2zl + £ a„z™, for 2 e B,(0) \ {0}, 



ia-i 

+ 

n=0 

with a_i € M, then the function 5ft(G) : K n -Bi(O) — > K has an analytic continuation on the whole 
disc B$(0). It follows from (|45[) that there exists a holomorphic function Q : B$(0) — > C such that 
Q(s) = 5R(g(s)) for any s € (-<5, J). Thus by Theorem 7.1, Chapter 4 of [25], we have 

3? (q( a )) e i7(«)et/2+v(-i7W)t ds = /" Q(s)e i7( S )et/2+v(-i 7 (,))t ds 



□ 



Letting e go to in equation (|42|) . applying Lemma [5.121 and the bounds developed above for the 
other integrals in (|42|) the theorem follows in the case x = 8/2. The case x = —8/2 is analogous. 
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APPENDIX 



A Proof of Lemma 12.61 

The proof is analogous to the proof of Theorem 12.21 The residue in Theorem 12.21 is equal to 1 (arising 
from the ^{-g/2<x<§/2} term), and from (|27p . the integral part is equal, for R large enough, reads 

(( [ +( ) Mz*L**d* \ , 

and the behaviour of these integrals follows exactly the lines of the proof of Theorem 12.21 



B Proof of Proposition 12.71 

Let us now consider a squared volatility of the form <j\ = a 2 + ai/t > 0, then the Black-Scholes call 
option reads 



n Cbs {S , S e xt , t, v t )=$l } ? Vt - e"* > Vt , (A-l) 



and let z± = (-x ± | (cr 2 + ai/t)) \ft/yja 2 + ai/i. Recall that (see [29] ) 

= 1 - §{z) = ° XP (~^ 2 ) (1 + (1A 2 )) , asz^+oo (A-2) 

ZV27T 

The case x > a 2 /2. As er 2 /2 = lim of/2, there exists <* such that for all t > t*, x > a 2 /2 = 
(a + ai/t) 2 /2. From (|A-1[) . we have, using a Taylor expansion for z±, 



(2^)- 1/2 cxp(-(VS s (.T,a) - x)t)A BS (x,a, ai ) (1 + 0(1/*)) 



The cases x < —a 2 /2 and — cr 2 /2 < x < er 2 /2 follow likewise. 
The case x = a 2 /2. From (|A-1I) . we have 

±», (*,*✓*«,«.*,) - « r-^=) - ^^±^ - <-^ts^ (1 + 0(1/i) ) 

<So V ' \vo- 2 i + ai/ art + (cHi + Oi)/2 

1 ■ "' /2 1 (1+0(1/0)4 + i2li^i(l + 0(l/ t) ). 



2 crV27rt crV27rt V27rf 
The case a; = — <t 2 /2 is analogous. 
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